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Abstract 

Many biological and artificial transport channels function without direct input of 
metabolic energy during a transport event and without structural rearrangements in- 
volving transitions from a 'closed' to an 'open' state. Nevertheless, such channels 
are able to maintain efficient and selective transport. It has been proposed that at- 
tractive interactions between the transported molecules and the channel can increase 
the transport efficiency and that the selectivity of such channels can be based on the 
strength of the interaction of the specifically transported molecules with the channel. 
Herein, we study the transport through narrow channels in a framework of a general 
kinetic theory, which naturally incorporates multi-particle occupancy of the channel 
and non-single-file transport. We study how the transport efficiency and the probabil- 
ity of translocation through the channel are affected by inter-particle interactions in 
the confined space inside the channel, and establish conditions for selective transport. 
We compare the predictions of the model with the available experimental data - and 
find good semi-quantitative agreement. Finally, we discuss applications of the theory 
to the design of artificial nano-molecular sieves. Key words: Transport; Channels; 
Selectivity; Efficiency; Diffusion; Occupancy 

1 Introduction 

The proper functioning of living cells involves continuous transport of various molecules 
into and out of the cell, as well as between different cell compartments. Such transport 
requires discrimination between different intra- and extra-cellular molecular signals 
and demands mechanisms for efficient, selective and specific transport (jlj). Specifi- 
cally, transport devices must be able to selectively transport only certain molecular 
species while effectively filtering others, even very similar ones. 

In certain cases, the selectivity and efficiency of the transport is achieved through 
direct input of metabolic energy during the transport event, in the form of the hy- 
drolysis of ATP or GTP ([l|). However, in many cases, molecular transport is efficient 
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and selective without the direct input of the mctabohc energy and without large scale 
structural rearrangements that involve transitions from a 'closed' to an 'open' state 
during the transport event. Examples of transport of this type include the selective 



permeability of porins 
plex in eukaryotic cells 




transport through the nuclear pore com- 

^ ^ artificial nano-channels and membranes, 

and other transport devices Ion channels iHI, [H, [H) , 
also belong to this class of transport devices - however the selectivity of ion channels 
depends on numerous factors that set them apart (13, 13) and place them beyond the 
scope of the present work. 

Transport devices of this type commonly contain a channel or a passageway through 
which the molecules translocate by facilitated diffusion. The selectivity and the ef- 
ficiency of transport are usually based not merely on the molecule size but on a 
combination of the size, strength of the interaction with the channel, speed of the 
spatial diffusion through the channel, and channel geometry (cf. Figs. 1 and 2). 



jatial diffusion through th^ 



H, [H, M, M, M, E M, M, Moreover 



large body of experimental data shows that the specifically transported molecules in 
many cases interact strongly with the channel (more strongly than the ones that are 
filtered out) and can transiently bind inside the channel (jJH, 0, H, 0, 0, B 0, U, 
[l8l.[l9l). Another important feature of such selective channels is that 
they are narrow, with a diameter comparable to the size of the transported molecules. 

Understanding mechanisms of the selectivity of transport through such channels 
is an important biological question and also has important applications in nano- 
technology and nano-medicine. For instance, it impacts creation of artificial molecular 
nano-filters. In addition, it poses a fundamental physical question: how does one make 
a selective channel that is always open, and does not have a movable 'shutter' specif- 
ically attuned to its corresponding transported molecules? Another important goal is 
to establish to what extent the theoretical models capture the essential properties of 
transport through narrow channels by comparing the models to experimental data. 

The precise mechanisms and the conditions for optimal selectivity of transport 
through such channels are still unknown. These systems span a wide spectrum of 
space and time scales and biological functions. For instance, porins are involved in 
the transport of small molecules into and out of the cell. They typically have channel 
length of several nanometers and a diameter of a couple of nanometers, tuned to the 
size of their corresponding transported molecules (e.g. water or small sugars) 
0,0,0)- The transport times through porins can be shorter than one millisecond |2)- 
In another example, the nuclear pore complex regulates the transport between the cell 
nucleus and the cytoplasm. It has a diameter of approximately 30 nm and a length of 

jass molecules up to 30 nm in size, within transport times of 
n. Artificial selective nano-channels have been constructed 



70 nm(Ii,i,0,E 
several milliseconds (155 
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5E 



several microns long and tens of nanometers in diameter that selectively transport 
molecules of sizes in the range of tens of nanometers. Such artificial devices have been 
used to selectively transport various molecules: including rnolecular enantiomers, short 
DNA segments and synthetic polymers dH, [13, H H O H El) • 



Nevertheless, it has been suggested that such channels might share common mech- 
anisms of selectivity and efficiency. Recent theoretical works propose a mechanism of 
selectivity that relies on two crucial factors, transient trapping of the cargoes inside 
the pore and the resulting confinement of the cargoes in the limited space within the 
channel. In particular , by modeling the transport as diffusion in an effective potential, 
the authors of (j3, HE Il3, S H [H, 113) have shown that the attractive interac- 
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tions of the transported molecules with the channel, such as transient binding of the 
molecules to binding moieties, increase the transport efficiency. More precisely, with- 
out an attractive potential inside the channel, the particles entering the channel have 
a low probability of traversing it to the other side. Attractive interactions inside the 
channel slow down the passage and increase the probability of individual molecules to 
translocate through the channel H, m, [13, [H, M, H Hi, [13, Ei). This mechanism 
of transport enhancement has also been known as 'facilitated diffusion' in the field of 
membrane transport 

However, space inside the channel is limited, and if the molecules spend too much 
time inside the channel, they prevent entrance of new ones. The channel thus becomes 
jammed and the tr ansp ort is diminished. To model the jamming of the channel, the 
authors of (|27l . [sol . l33l . [s^, Issh assumed that additional molecules cannot enter the 
channel already when one molecule is present inside. They showed that particles 
whose interaction with the channel is weaker than the optimal, have a low probability 
of traversing the channel, while particles that interact too strongly with the channel 
jam the transport. This allows discrimination between the molecules based on the 
strength of their interaction with the channel, and provides a mechanism of selective 
transport; transmission efficiency is optimized for a particular interaction strength and 
rate of transport. Optimal trapping time, which maximizes the transmitted current, 
has also been demonstrated for single-file transport in (j48l ). 

However, during transport, the channel can be occupied by multiple molecules, 
which cannot bypass each other, or do so only in the limited fashion, due to the 



confinement in the limited space inside the channel y, |45|, |47|, |48|, S, [53,[5ll). The 
transport is not necessarily single-file: the number of molecules that can be present 
at a position along the channel depends on the ratio of the channel diameter to the 
molecule size. We must also recognize that the transport properties of narrow channels 
are not dominated by the equilibrium thermodynamic channel-cargo interactions per 
se, but by the rates at which the cargoes enter, translocate through and exit from 
the channel with a potentially complicated geometry (Q, S, i S m, 0, M, 0, Si). 
For instance, the trapping time in the channel can be limited by the time it takes 
to find a narrow exit from the channel by diffusion. This phenomenon is known as 
entropic trapping In the case when the rates are determined solely by 

the interaction strength, stronger interactions with the channel imply slower rates and 
higher trapping times (cf. Figs. 1 and 2)(|i, S S, [H, 53). 

Understanding the effects of multiple channel occupancy and jamming on the trans- 
port selecti vity is especially pertinent to the analysis of single molecule trac king ex- 
periments ([57, H, 53) and the design of artificial nano- molecular filters ©Hill, 



In this paper, we analyze transport through narrow channels in a framework of 
a general kinetic model based on exclusion process theory as a function of the ki- 
netic parameters of transport Specifically, we examine the rates of entrance, hopping 
through and exit from the channel. We extend the previous work to include multiple 
occupancy and inter-particle interactions inside the channel beyond single file. We 
investigate how the concentration of the cargoes, the channel length and radius, the 
dimensions of the transported molecules and the interactions between them inside the 
channel influence the transport. An important goal of this paper is to explore whether 
a theory that has only two essential ingredients: 1) transient trapping of the molecules 
inside the channel and 2) inter-particle crowding due to the confinement in the lim- 
ited space inside the channel, can provide an adequate explanation of the selective 
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transport through narrow channels by comparing the predictions of the theory with 
the available experimental data. 

The paper is organized as follows. We first discuss a channel that consists only 
of one 'site' and then two-sites. Next, we discuss transport in a uniform symmetric 
channel of arbitrary length, for both single- file and non-single- file transport, and es- 
tablish conditions for optimal transport. We then discuss the transition between two 
transport regimes, jammed and un-jammed, and establish the relative contribution of 
the jamming of the channel entrance as compared to crowding inside the channel, to 
the transport selectivity and efficiency. Next, we compare predictions of the theory 
with the experimental data. We conclude with discussion of the results, their relation 
to the previous work, and consider potential applications. 



2 Inter-particle interactions inside the channel 



A transport channel can be regresented as a chain of positions (sites), as illustrated 



in Figs. 1 and 2 |27|, l29|, |35|, \M \M, \M |4g|, l50|, |51D- The particles attempt to 



enter the channel at a given position, with an average rate J and subsequently hop 
back and forth between adjacent sites, if those are not fully occupied, until they 
either reach the rightmost or leftmost sites, from where they can hop out of the 
channel. Hopping out from the rightmost site represents the particle reaching its 
destination compartment, while hopping out from the leftmost site channel represents 
an abortive transport event, where the molecule does not reach its destination (cf. 
Figs. 1 and 2). In the continuum limit, when the distance between the adjacent 
sites tends towards zero (and their number to infinity), with an appropriate choice 
of the transition rates between the sites, the problem can be reduced to diffusion in 

Si 13, |43 ) (cf. also Appendix) . Note that the 



an effective continuous potential ()27l . [2; 



discrete positions (sites) do not represent the actual binding sites inside the channel. 
Rather, they are a convenient computational tool that allows one to explicitly take 
into account competition for space and interactions between multiple particles inside 
the channel {2^ 2^, 13, 48, 4^, 5^, 51)- The distance between the positions reflects 
the size of the particles. 

As the particles accumulate in the limited space inside the channel, they start to 
interfere with the movement of the neighboring particles and prevent the entrance of 
new ones. We must differentiate between the speed, the efficiency, and the probability 
of transport. The speed is determined by the time the particles spend in the channel. 
The efficiency of transport is determined by the fraction of the impinging flux that 
reaches the rightmost end. It depends on the kinetic parameters of the channel, such 
as transition rates inside the channel and the exit rates at its ends. The selectivity of 
transport is determined by the different efficiencies at different values of the kinetic 
parameters (26, 27, [2^, [sO, 31, 33, [ish . Transport efficiency is different from the 



probability that an individual particle translocates through the channel. The latter is 
defined as the fraction of the particles that reach the exit after entering the channel. 
We discuss these issues in detail below. 



2.1 One site channel 



To get started, let us consider a 'one-site' channel (|27l . l42l ). where all the internal 
spatial and energetic structure of the channel is absorbed into the forward and the 
backward exit rates r^, r^. 
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Kinetic diagram of such a 'one-site' channel is shown in Fig. IB. The state of the 
channel is specified by the particle density (0 < n < 1) at the channel site (or, in other 
words, the probability of the channel to be occupied). It obeys the following kinetic 



equation (|27l.l42f): 

h = J(l — n) — (r^ + r_+)n (1) 

which takes into account that the particles can enter the channel only if it is not 
occupied. The average time a particle spends inside the channel is r = l/{r^ + r^) 



At steady state [h = 0) we get for the average density and the forward flux: 

(2) 



Jo 



J + J+ 1/r 



As mentioned above, we define the transport efficiency as the ratio of the transmitted 
flux to the entering flux, Eff = Jout/J- Thus, from the equation ([2|) we learn that 
the transport efficiency Eff(J, r^,r^) = — is a monotonic function of both 

the forward exit rate r_> and the backward exit rate r^. Therefore, for the 'one-site 
channel' there is no optimal combination of the exit rates that would maximize the 
transport. As we shall see, this is not the case for longer channels. However, even a 
single site channel can have more interesting behavior, if the forward and the backward 



exit rates are not independent (|27l . l33f ) (cf. Appendix) 



2.2 Two site channel 

Let us consider now a longer channel consisting of two sites: 1 and 2. This is the 
shortest channel that explicitly takes into account the asymmetry between the channel 



entrance and exit, and exhibits non-trivial transport properties ()4l.l27ll33l.l34l.l49l). The 
kinetics of transport through such a channel is illustrated in Fig. 1 C and D. The 
particles enter the channel at the entrance site 1 with an average flux J, if it is un- 
occupied. The backward exit rate to the left from site 1 is and the forward exit rate 
to the right from site 2, is r^. Once inside, a particle can hop back and forth between 
sites 1 and 2 with rates ri2 and r2i, respectively, if the target site is un-occupied. The 
exit rates and can be thought of as the 'off' rates for the release of the particles 



from the channel (j27T ). For simplicity, in this section we assume that the channel is 
internally uniform and symmetric with ri2 = r2i = r and = r<_ — and that 
each site can be occupied only by one particle. 

The state of the channel is characterized by the average occupancies of the sites, 
< ni < 1 and < 712 < 1- For an internally uniform channel, these average 
occupancies can also be viewed as the probabilities that the sites 1 and 2 are occupied 
by a particle (0, S, I47I ). The kinetic equations describing transport through such a 
channel are (Fig. IC) 

rii = J(l — Til) - ToTii + r7i2(l — ni) — rni(l - n2) (3) 
«2 = rni(l - 712) - ^2(1 - rii) - ron2 

and the transmitted flux is Jout — 7'o'^2- 

The transport efficiency Eff(ro) = Jout/J is the fraction of the flux J that exits 
the channel to the right. Solving equations Q at steady state (rii ~ 712 = 0), one gets 
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for the transport efficiency: 

TV 

Importantly, unlike in the one-site case, for a given entrance flux J, the transport 
efficiency Eff(ro) has a maximum at a certain value of the exit rate r"^^^ = ^/Jr. This 
provides a mechanism of selectivity; only particles whose residence time in the channel 
(determined by the interactions of the particles with the channel) is close to 1/ j-™^^ 
are transmitted cfliciently (jl^, 0, 0, [sj, [H, S) . 



The total efficiency Eff = Jout/J is influenced by two different effects, the jam- 
ming of the channel entrance and the mutual interference between the particles inside 
the channel. The flux that actually enters the channel is Jin = J(l — n,i). The re- 
maining portion of the flux, Jni, does not enter the channel because the entrance 
site 1 is occupied ni fraction of the time. The fraction of the entering current Jin 
that reaches the exit on the right determines the transport probability = Jout/Jin, 
which characterizes transport through the channel. From the equations ([3]) 

T 

= TT-— (5) 

2r + To 

Very importantly, is independent of the flux J and is equal to the efficiency in 
the single-particle limit, J ^ 0. That is, it is equivalent to the probability of a single 
particle to translocate through the channel when no other particles are present. Thus, 
surprisingly, the crowding of the particles inside the channel does not, on average, 
influence their movement through the channel. We discuss this effect at length below. 

To summarize this section, selective transport can arise from a balance between 
two competing effects, enhancement of the transport by the transient t rapping and the 



eventual jamming of the channel if the trapping times arc too high ()27l .l2 



appmg anc 



2.3 Channel of arbitrary number of sites 

In this section we study transport through a channel of arbitrary length, which is 
modeled as a chain of N positions (sites): 1, 2...i...N. Particles enter at site M (not 
necessarily the leftmost one) with an average flux J if the entrance site is not fully 
occupied. Once inside the channel, a particle at site i can hop to an adjacent site z± 1 
if the latter is not fully occupied. In order to model the excluded volume interaction 
between the particles in the channel (they can bypass each other only in a limited 
fashion) , we introduce the maximal site occupancy n,„ , which depends on the ratio of 
the channel diameter to the size of the particles. When at an outermost site 1 or N, 
a particle can leave the channel with the rate and r_,, respectively, or hop into 
the channel with the average rate ri^2 or tat^jv^i, respectively. The kinetics of this 
process is illustrated in Fig. 2 A. The rates ri^i±i determine the speed with which 
the particles diffuse through the channel, while the exit rates and r_, reflect how 
fast the particles can leave the channel. The kinetic rates ri^i±i, r^ and r^ are 
determined by the microscopic interactions of the particles with the channel, and by 
its geometry. As before, the exit rates r_> and can be thought of as the 'off' rates 
for the release of the particles from the channel ((27h . In general, with a proper choice 
of the transition rates, ri^i±i / ri_^i±i = exp(t/i_|_i — Ui)j2) in the continuum limit 



the model reduces to diffusion in the potential U{x) (j27l . |32|. |41|. l4 2| ). Tra ppin g of the 



particles in the channel corresponds to low exit rates r^, < r (27, H, 4^, 4 
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For simplicity, we assume that the channel is internally uniform, such that all 
the internal transition rates are equal, ri_,i±i — r for all i. At any time t, the 
state of the channel is specified by the number densities of the particles at each site 
TT-i, 712, n^, ...njv. The kinetics of transport through such a channel is described by 
the following equations (Q, H [i^, [H, Si, [HO) 



) + mi„i(l ) + rui 

Tim. 



l(1 



rni{l 



n,- 



rni{l 



= JSi,m{^ ) + r{ni-i + ni+i - 2ni) 

with the boundary conditions at sites 1 and A'^ 



rii = J(5i,m(1 - — ) 

Tlnr 



-ni 



rni(l ) 

Tim, 



rn2{l ) 

Tljn. 



(7) 



UN 



n. 



UN 



™jv-i(l ) = -{r + r^)nN + rriM- 



J6i.m{1 ) - (r + r^)ni + rn2 

—r^ntq — rn]y{l — 

where the (S-function is Si_j = 1 if i ^ j and zero otherwise. The terms ni(l — 
ni±i/n„i) in equations ^ and ([7]) reflect the fact that a particle can jump to the 
next site only if it is not fully occupied, 7ii±i < n„i. Importantly, for an internally 
uniform channel, at all the internal sites the cross-terms of the form nini±i cancel out 



(471. ksl). For such uniform channels, the equations ^ and ([71) are exact and ni/rim 
is equivalent to the probability of a site i to be occupied (|47l . |48I ) . Obstruction of the 
space inside the channel by the particles inside of it affects only the entrance to the 
channel at site M . 

We define the efficiency of transport as the ratio of the forward exit current Jout — 
r—tJiN to the incoming flux J, Eff(ro) = Jout/ J- It is the fraction of the incoming flux 
that traverses the channel. Note that the efficiency is different from the probability 
of individual particles to traverse the channel after they have entered, because some 
of the particles attempt to enter the channel and are rejected if the entrance site is 
occupied. 

The linear equations (H]),® can be solved analytically for any N (fish . In the fully 
symmetric case, when the forward and the backward exit rates from sites 1 and N are 
equal, = = Tq, the efficiency is given by: 



Eff(ro) = 



(r + (M - l)ro)r„ 



ro(2r +{N- l)r„) + ^{r + {N - l)ro + (M - 1){N - M)rl/r) 



(8) 



(for M < N/2). 

Note that in the single particle diffusion limit, J ^ 0, the efficiency Eff M/ {N + 
1) without trappin g (r- p —>■ r) and Eff 1/2 for strong trapping {rg ^ 0), in accord 
with known results 13, [13, SI [41I. |4^. [isl. [55I) . Essentially, without transient trapping, 
the probability to traverse the channel is low and reaches one half for very strong 
trapping. 

One may rewrite equation ([8]) in terms of the trapping time, which is equal to 
r = 27- iHSliS) to arrive at 



Eff(r„) = 



(2rr + (M - l)N)N 



{N - l)N + Arr + ^{jjirr? + {N - l)Tr + [M - 1){N - M)) 



(9) 



H+l ■ 



(6) 
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where J' = J/r is the normaUzed flux. Note that the transport efficiency does not 
depend on the absolute values of the transport rates r and Tq, but only on the nor- 
malized parameters rr and J/r. This means that the transport efficiency can be the 
same for different particles, even if the kinetics of their transport through the channel 
is very different from each other, as long as they possess the same rr and J/r. 

As already seen in the two-site case, the transport efficiency Eff (ro) of equation 
([8]) has a maximum at a certain value of the exit rate (for M = 1) of 



rrVr^^^§^ (10) 
and the maximal flux at this rate is 

((TV- l)r^a^/r + l)Vl 

(cf. Appendix for M ^ 1). 

This feature provides a mechanism of selectivity; only the particles whose exit rate 
is close to the optimal one, r™^'^, are transmitted efficiently. Particles with exit rates 
higher than the optimal have a higher chance of returning back because they don't 
spend enough time inside the channel in order to reach the farther exit on the right 
side. On the other hand, due to the limited space inside the channel, the particles 
with the exit rates lower than optimal spend so much time in the channel that it gets 
jammed and the entrance of new particles is inhibited 

Equations ([5]) and pT7)l qualitatively agree with the results of (30|, |33|, |3J, l35r ). 
which assumed that only one molecule can occupy the channel. Figure 3 shows how 
the transport depends on the channel length N, the entrance flux J, the exit rate 
and the effective channel width n„i. Note that the optimal exit rate decreases with 
the channel length N; for longer channels, a particle has to spend more time in the 
channel in order to reach the other end. Also note that the optimal rate of equation 
(|10p . r^^^/r, is less than one for J/r < TV — 1; the optimal interaction is attractive 
for small currents and long channels. We elaborate on this issue in Appendix. 

2.4 Transport efficiency vs. translocation probability 

In this section, we elaborate on why the flux through the channel decreases in the 
limit of very low exit rates? Is this because new particles cannot enter or because the 
particles inside the channel interfere with each other's passage? 

The fraction of the incoming flux J that actually enters the channel is Jin = J(l — 
ni/rim). The remaining portion of the flux J^ cannot enter because the entrance 
site is occupied on average ^ fraction of the time (cf. Sec. 4 for calculation of the 
densities). The total efficiency is determined by two quantities: i) the fraction of the 
flux that enters the channel Jin and ii) the fraction of the particles that upon entering 
the channel, actually reach the rightmost end. The latter defines the probability = 
Jout / Jin of a particle exiting to the right after it has entered the channel and is given 

by 

P ^ -^""t ^ r + (A/ - l)ro 
^ Jin 2r + {N-l)r, ^ ^ 

Remarkably, it is independent of the flux J and is exactly equal to the efficiency in 
the single particle transport limit, J — > 0. This means that in unform channels the 
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interactions between the particles in the channel do not affect the transport probabil- 
ities of the individual particles. The effect of the channel occupancy manifests only in 
the jamming at the entrance site. 



3 Optimal transport and jamming 

The lower the exit rate Tq, the longer the time that the particles spend inside the 
channel. The trapping time varies as r = ^7- (|27. . 43, 52). As shown in the previous 
section, at very small exit rates Tq, the trapping time is so high that the channel 
becomes jammed. Thus, the transport efficiency is maximized at the particular exit 
rate r'^^^. Inspection of the Fig. 3A reveals two distinct transport regimes, roughly 
separated by the maximum of the transport efficiency at Tq = r^^^. At the high 
values of Tq > r™'*'^ the transport of individual particles is essentially unhindered 
by the presence of the others, as evidenced by the fact that the transport efhciency 
curves collapse onto the dashed line, representing the zero-current, single-particle limit 
(Fig. 3A). At the low values of the exit rate where Tq < r™'*'^, the accumulating 
particles start to obstruct the entrance of the new ones. This feature provides a 
natural definition for the 'jamming transition' around the = r™^'^. 

Solving equations ^ and (O, we get for the density profile of the particles inside 
the channel, at the steady state: 

J{r+{N-i)ro) ^^g^ 



roi2r+{N^l)ro) + Jir + {N-l)ro) 

(for il/ = 1, n,„ = 1). Note that unlike the equilibrium distribution, the maximum of 
the density profile is near the channel entrance at site 1. 
The total number of the particles in the channel is 

N -V^ _jV J{2r + (N - l)r,) 

iVtot - - 2 ro(2r + {N - l)r„) + J(r + (A^ - l)ro) 

Note that in the limit j-q 0, A^tot ^ A^- That is, the particles accumulate and never 
leave the channel. Therefore, from equation (fTil) one finds that at the point of the 
jamming transition, = r^^^^ the number of the particles in the channel is 



^jam = N— -X^i (15) 

2(j/r + 2(^+^)j 

Equation ([T5|) has important consequences (cf. Fig. 4 for illustration). It shows that 
for long channels, where A^ 3> J/r, the fraction of the occupied sites at the jamming 
transition tends to one half: Ajam/A 1/2 when A^ ^ J/r. This means that 
long channels can be filled up almost to half of their maximal capacity A^ before the 
jamming effects start to matter. For the occupancies below the jamming transition, 
the particles travel through the channel essentially unhindered. These effects are 
illustrated in Fig. 3A and Fig. 4A. This might explain why experiments on transport 
through narrow channels often measure apparent diffusion coefficients that are almost 



as large as those for the free diffusion (|53l . |5E 
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3.1 Jamming and saturation of the flux through the channel 

Although the transport efficiency Eff (ro, J) decreases with the increasing flux J, the 
total transmitted flux Jout = J EfF(ro, J) saturates at large fluxes (J/r oo) to the 
limiting value 

(for M = 1). This saturation of the transmitted flux at large incoming flux J is 
another manifestation of the jamming of the channel entrance by the particles inside. 
Indeed, equation ([T^ shows that the density at the entrance ni tends to tii = 1, as 
J/r — > oo. In other words, the flux saturates because no more particles can enter the 
channel. This is neatly summarized by the observation that ni = Jout/J^t- 

By contrast, the exit site N is not completely blocked even at high J and npf 
1/(1 + {N — l)ro/r) as J — > oo. Thus, even at very large fluxes, when the entrance 
site is completely blocked, the channel is not fully occupied. From equation (|15p . the 
number of particles in the channel is 

^ N2r+{N^l)r^ 
'°* 2 r + {N- l)ro 

In particular, for long channels (A^ — 1 3> r/ro), the channel occupancy in the saturated 
limit is Ntot/N = 1/2. Also note that the saturated flux is proportional to Um, and 
that it decreases with rg/r. 

The results of this section closely parallel Michaelis-Maenten kinetics of multi-step 
enzymatic reactions (|46l ) and are important for the estimation of binding affinities 
from the channel transport experiments (5^, 54, 5^, as well as for comparison with 



experiments on flux through artiflcial nano-channels ( - Sec 4). [3) 



4 Comparison with experiments 

In experimental systems, the exit rates and the rates of transport through the channel 
are determined by a potentially complicated kinetics of binding and unbinding in- 
side the channel. Can the theory adequately describe these experiments? Facilitated 
diffusion theories produced results consistent with the experimental observations of 
the transport of gases through functionalized membranes [stI . [39! ) enhancement of 



transport of oxygen by myoglobin (|38l ) and the transport through bacterial porins(jj). 
In this section, we compare the theoretical predictions of this paper with the experi- 
ments of Ref. (HI). 

Briefly, in the experiments of Ref.lflih that we chose for comparison with theo- 
retical predictions, transport of short DNA segments through artificial nano-channels 
was studied. The flux of the DNA segments through these channels was measured in 
two cases: 1) empty channels and 2) channels were lined with single-stranded DNA 
hairpins, grafted to the walls. Each hairpin has a stretch of 18 un-paired bases in 
the middle. The transported particles were 18 base ssDNA segments with the se- 
quence complementary to the un-paired regions of the ssDNA hairpins inside the 
channels. Thus, the transported DNA segments can transiently hybridize with the 
DNA grafted inside the channel. That investigation found that the flux through the 
DNA-containing channels is higher than through the channels without DNA hair- 
pins inside, providing evidence that the transient trapping indeed facilitates transport 
through nano-channels. However, eventually the interactions between the particles in 
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the limited space inside the nano-channel block the passage, causing the transmitted 
flux to saturate with the increase in the incoming flux. This is another signature of 
the transient trapping discussed above in Sec. 3 flQElii). 



The radius of an empty channel is i? ~ 6 nm and the channel length is 
The grafted ssDNA hairpins reduce the passageway radius, which, for the purposes 
of comparison with the theory, we roughly estimate as i? ~ 3 nm for the channels 
with DNA grafted inside. Using the value for the gyration radius of the transported 
DNA segments 5^1 nm 

(0, 

md ). we estimate n,„ = 6 for the empty channels and 
rim = 3 for the channels with the grafted DNA hairpins inside. Furthermore, we 
estimate the incoming flux as 4Z)outci? (41), where I?out — fitnJ,, diffusion 
coefficient of the transported DNA coils outside the channel md): 77 is the viscosity of 
the solvent, Sh — 0.75 is the liydrodynamic radius of the coils (|60l) and c is the outside 
concentration of the transported DNA (|4ll . |42| ). To model the finite capacity of the 
channel, we estimate the number of available positions in the channel as iV = L/{2S), 
where L = 6^ (0, S, 0, [5ll)(cf. also Appendix). 

Finally, Vo/r = ^^^j^j^Z, where Z is the reduction in the exit rate due to the 
trapping inside the channel (|27l )(cf. also Appendix). We return to the question of 
how Z is related to the actual binding energy below. The ratio Do^t / ^in and Z are 
the two independent fitting parameters of the model (note that the r and Tq appear 
as independent parameters in eq.®. 

We first tested the model for the case without DNA segments attached inside 
the channel. The data (black dots) and the fit (black line) with £'in/£'out = 0.42 and 
Z ~ 1 are shown in black dots in Fig. 5A. Analogously, for the channels with the DNA 
hairpins inside, the fit of the equation ([8]) to the data (red dots) is shown in the red 
fine in Fig. 5A, with the best fitting parameters An/-Dout = 0.004 and Z = 0.00007. 
Note that the diffusion coefficient is significantly reduced in the narrow channel filled 
with the grafted DNA hairpins ([H). 

As expected, the transient binding of the transported DNA segments to the DNA 
hairpins inside reduces the exit rate Tq by a factor Z = 0.00007. Because this energy 
is influenced by many factors that are poorly understood 

13), 

one can not easily 

connect the value of Z to the actual binding energy between the transported DNA 
and DNA hairpins. However, if the reduction in the exit rate is indeed determined 
mainly by the effective binding energy e, then the function Z ~ exp(— e/fcT) should 
describe the trend in the dependence of Z on e 27, [1^, \^ [ish . The authors 
of (flih measured fluxes through the channel for DNA segments possessing different 
numbers of mismatches to the DNA grafted inside and found that the flux decreases 
with the number of mismatches. Thus, assuming as a first approximation that the 
binding energy e decreases linearly with the number of mismatches n, so that e{n) ~ 
e„=o(18 - n)/18 we get Z = exp(ln(0.00007)(n - 18)/18). The prediction of equation 
((8)) with this choice of Z is compared with the data in Fig. 5 B. It shows that this 
simple estimate correctly reproduces the trend in reduction of the transmitted flux 
with the number of mismatches. Note that there are no additional fltting parameters 
used in this figure. 

That the simplified theory developed in this paper can correctly reproduce the 
trends in the observed fluxes, and even gives semi-quantitative fit of the data for 
reasonable values of the parameters, is encouraging. This demonstrates that a theory 
that is built upon only two essential assumptions, 1) facilitation of diffusion by the 
transient trapping inside the channel and 2) mutual interference between the particles 
crowded in the confined space inside the channel, does provide an adequate explanation 



Selective transport in nano-channels 



of the experimental data. Moreover, the theory provides verifiable predictions about 
how the flux should change with the channel diameter and length, as well as the particle 
size and concentration, as described in section 3. Comparison of these theoretical 
predictions with future quantitative experiments will lead to further ramification of the 
theoretical approach and will facilitate the design of artificial selective nano-channels 
with desired properties. 

Discussion of other effects observed in (|l6t ). which are attributable to a confor- 
mational transition of the hairpin layer during transport, is outside the scope of the 
present work. 



5 Summary and discussion 



Proper functioning of living cells requires constant transport of different molecular 
signals into and out of the cell, as well as between different cell compartments. To 
carry out this task, the living cells have evolved various mechanisms for efficient and 
selective transport. 

One class of transport devices comprises narrow channels whose diameter is com- 
parable to the size of the molecules transported through it. Examples include selective 
transport through the nuclear pore complex, bacterial porins. 111,!!, 0, 3, [ p . 0J3, H 
Tol U, 12) and other non-biological transport systems such as zeolites P5ri40ll49l). A 
crucial feature of such transport channels is their ability to selectively transport their 
specific signalling molecules while efficiently blocking the passage of all others. 

Driven by the notion that natural evolution has optimized the function of such 
devices, large effort is being currently invested into the creation of artificial nano- 
molecular sorting devices that mimic the function of biological channels 111, Q M, 
3 E M, H, [Hi). The desi gn of such devices requires detailed understanding of the 



principles of selective transport through narrow channels. 

The precise conditions for the optimal transport selectivity through narrow chan- 
nels still elude our understanding. A large body of experimental work indicates 
that the selectivity is often based on the differential interactions of the transported 
molecules with their corresponding transport channels. Moreover, interaction of the 
transported molecules with their corresponding transport channels is strong, exceed- 
ing the interaction of the non-specific competitors. Another salient feature of such 
channels is that they are narrow, so that the particles cannot freely bypass each other 

(2, 1 i, i, i S i i M M M 



Recent theoretical works have shown that selective transport through narrow chan- 
nels can arise from a balance between efficiency and speed; transient trapping inside 
the channel increases the probability of a molecule topass through the channel, but 
leads to jamming at too high trapping times li,!!!,^,^^,!!!!,!!!,^,!!!,^,!!^,® 

Extending previous work, in this paper we have analyzed transport through narrow 
channels in the framework of generalized kinetic theory. We represent a transport 
channel as a sequence of positions (sites) and the transport through the channel is 
determined by the hopping rates from one position to another inside the channel, as 
well as by the exit ('off') rates from the channel at its ends. To take into account 
the limited space inside the channel, and the finite size of the transported particles, 
we allow only a limited occupancy at each position, Um- Thus, a particle present 
at a given position along the channel, can hop to an adjacent position only if the 
latter is occupied by less than rim particles. Our model allows one to naturally treat 
channel occupancy by multiple particles and extends the treatment beyond the single- 
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file transport. The main determinant of tlie transport properties of the channel is 
not the interaction strength of the particles with the channel per se, but the kinetic 
properties of the channel, which determine the trapping time r and also depend on 
the geometrical properties of diffusion in the confined space inside the channel. These 
possibilities are illustrated in the Fig. 2. 

We briefly summarize our major findings below. In qualitative agreement with 
previous work, we find that the transient trapping of the particles in the channel 
increases the transport probability; particles that have high exit rates do not stay in the 
channel long enough to reach the exit into the destination compartment on the right 
side (cf. Fig. 2), and have a high probability to return back HI, [11, [3l|, [13, Si) . 
Essentially, transient trapping increases the time that the particles spend inside the 
channel to be long enough in order to reach the exit on the right side. Thus, although 
each individual particle spends more time in the channel, the transmitted flux is 
higher. If the only measured quantity is the flux through the channel, experiments 
can not easily distinguish between the probability of transport and transport speed. 
H, [Hi, [Hi, [Hi). However, they can be distinguished in the experiments that follow 
transport of individual molecules ([HtI . [HsI . [Hoh . 

When the exit rate is too slow or the incoming flux is too high, the rate of particles' 
entrance to the channel becomes higher than the rate of exit and the particles start 
to accumulate inside the channel, because the space inside is limited. This leads to 
two distinct effects. First, the particles inside the channel start to interfere with the 
passage of each other. Second, they block the entrance site and inhibit the entrance 
of new particles. The channel thus becomes jammed. We must distinguish between 
translocation probability and transport efficiency. Transport efficiency is the fraction 
of the total incoming flux that reaches the exit. Only a certain fraction of the in- 
coming flux can enter the channel because the entrance site can be occupied when 
particles attempt to enter. This effect decreases the capability to enter the channel 
and, as a consequence, decreases the transport efficiency. Interactions between the 
particles inside the channel can also influence the probability of individual particles to 
translocate through the channel upon entering compared to the single particle case. 
However, we found that for internally uniform channels the crowding of the particles 
inside the channel does not affect the probability of individual particles to translocate 
through the pore. Thus, the effect of particle accumulation in the channel manifests 
only in the blocking of the entrance to the channel, which leads to the decrease in the 
total transport efficiency (and transmitted flux) at low exit rate or high incoming flux 
(Fig 4B). Thus, we predict that the kinetic profile near the entrance is an important 
factor in determining the selectivity of transport. 

For symmetric channels, this balance between the transport probability and the 
obstruction of the particle entrance to the channel, determines the optimal exit rate 
^max^ (cf. Fig. 4) which maximizes the transport. This provides a basis for selectivity, 
whereby different molecules can be selected by the kinetics of their transport through 

the channel ([23, [H, [13, [Hi, [H, [Hi, [ii) . In the case discussed in this paper when many 



particles can be present in the channel simultaneously, the optimal exit rate and the 
optimal flux depend on the length of the channel (cf. Section 3). Notably, this is a 
purely kinetic selectivity mechanism: although a low exit rate can be due to energetic 
interactions between the transported particles and the channel; the transport efficiency 
is not determined by the equilibrium occupancy considerations; the selectivity can go 
beyond the difference in the equilibrium binding affinities between different molecules. 
The fact that the transport efficiency has a maximum at a certain value of the exit 



Selective transport in nano-channels 



rate = y-max pj-Qvides a natural definition for the 'jamming transition'. Particles 
with the exit rates faster than r™*^^ pass through the channel essentially unhindered 
by the interactions with other particles because they do not stay in the channel long 
enough(Fig. 3A.) On the other hand, particles with exit rates slower than rmax 
compete with each other for entrance into the limited space inside the channel and 
the channel becomes jammed. Importantly, we found that the interactions between 
the particles, and the competition for the limited space inside the channel do not 
play an important role until quite a few of them accumulate in the channel. For 
long channels, approximately half of the available channel sites are occupied at the 
jamming transition ( Fig. 4). This implies that in many experimental situations the 
interactions between the transported particles do not play a significant role, and may 
explain why the apparent diffusion coefficient in many flux measurement experiments 
is found to be almost as high as for free diffusion (js, 53, 55, 58). 

Although many particles can be crowded inside the channel, and the entrance to 
the channel is blocked, transmitted flux does not disappear even at high fluxes and 
densities, but rather saturates to the limiting value determined by the trapping time 
and the channel length (cf. Fig. 3 and Fig. 5.) This closely parallels Michaelis- 
Maenten kinetics of enzymatic reactions ([l|, l46h and mig ht be relevant to estimation 
of binding strengths from flux experiments (jSSl. [sl . l56l ). Notably,the selective and 
efficient transport persists beyond the single file transport, even when the ratio of the 
channel diameter to the particle size is large. In this case, the optimal exit rate r™'^^ 
is simply shifted to lower values. 

In order to determine whether the theory developed in this paper can provide 
an adequate description of experiments, we compared predictions of the theory to 
the experiments reported in (|l6l). That work found that at low concentrations of 
the transported particles the flux through artificial nano-channels increases if the 
particles can transiently bind inside the channel. Moreover, as the binding energy of 
the particles was decreased, the enhancement of the flux was lower. However, as the 
concentration of the particles in the origin compartment increases, the flux saturates 
for the channels with transient binding, while the saturation if not observed for non- 
binding channels, at the experimental range of concentrations. Both these results 
are in agreement with the theory and can be semi-quantitatively described by the 
theoretical predictions, as shown in section |4l 

Thus, wc find that the theory based on only two main ingredients: 1) transient 
trapping of the molecules inside the channel and 2) crowding of the molecules in 
the limited space inside the channel, captures the essential features of the selective 
transport through nano-channels. Moreover, the theory provides verifiable predictions 
regarding how the flux and selectivity of such channels depend on the channel length, 
channel radius, the size of the transported molecules and the strength of the interac- 
tions of the molecules with the channel. In particular, we predict that the flux through 
such channels can be optimized by varying the interaction parameters and the chan- 
nel dimensions. Such predictions are useful for the design of artificial nano-sorting 
devices. Further quantitative experiments and comparison with the theory are needed 
in order to test the theory and for its further refinement. 

We expect that the effects described in this paper should play a role in selective 
transport through any narrow channel. For instance, the effects described in this 
paper might be relevant in determining the selectivity of the ion channels, although 



other factors might be dominant (|2ll . |22| . |23| . 124 ). In each particular system other 



effects related to molecular details might be dominant determinants of selectivity. 
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Such effects might include the long range electrostatics and channel fluctuations in 
the ion channels, the details of the transfer of the transported molecules from one 
binding moiety to another, and conformational changes of the filaments that carry the 
binding moieties (as in the nuclear pore complex and other polymer-based systems). 

Finally, we note that the theory developed in this paper can also be applied to 
other signal-transducing schemes, such as signalling cascades and multi-step enzymatic 
reactions (jH, 31, .62, 63). 

The author is thankful to C. Connaughton, B. Chait, I. Nemenman, J. Pearson, A. 
Perelson, Y. Rabin, K. Rasmussen, M. Rout, N. Sinitsyn, T. Talisman, Z. Schuss for 
stimulating discussions, P. Welch for comments on the manuscript, and anonymous 
reviewers for helpful suggestions. This research was performed under the auspices of 
the U.S. Department of Energy under contract DE-AC52-06NA25396. 



Appendix 

Single particle occupancy: connection to previous work 

In this section we show that the model of this paper can be reduced to previous 
models, in a proper limit. Let us assume, following (27, 13, 33, 13, 35) that already 



when the channel is occupied only by one particle, it prevents the entrance of others. 
The channel, however, is long, and the particle can obey complicated kinetics inside, 
which determines its probability to traverse the channel, and the time it spends inside. 
Physically, such situation can arise, for instance, due to strong long-range repulsion 
between the particles. 

In this case, the problem reduces to a 'single-site' channel of Sec. 12.11 but with 
forward exit rate r_^, backward exit rate that are not independent, but are de- 
termined by the internal kinetics of the channel, and are related through the single- 
particle dwelling time r and transport probability Ptr- As in Sec. 12. 1[ the transmitted 
flux is 

Jout = -rr^ — (17) 

From equation ([5]), the probability of a single particle to traverse the channel of 
length TV (for J ^ 0) is Ptr = 1/(2 + (A^ - l]iV/(rr)) = r^/{r^ + r^), and the 
residence time is r = N/(2ro) = ^/{r^ + r^) (|52l ). Thus, we get for the transmitted 
flux: ^ 

ni + jr) (l + ^) 



which is identical to expressions obtained in Ref.(j27(), if one bears in mind that the 



flux is J = fconC, where c is the concentration of the particles outside the channel. 



In is important to note that the optimal exit rate in this case is r^^^^ = ^ , 
that is almost independent of A^ for long channels. This is in contrast to the model 
of Sec. 2.4, which takes into account multiple occupancy of the channel by many 
particles - where the optimal exit rate decreases with N . The optimal current is, by 
contrast, higher for niultiple-occupance channels. This is natural - if more particles are 
can occupy the channel before it becomes jammed, the channel can sustain a higher 
current. 
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Connection between continuum and discrete models. 



Discrete model of equation ([6]) reduces to a continuum description of transport inside 
the channel, if one defines the one-dimensional particle density c^{x) = rii/a where 
a is the distance between the 'sites', so that x = ai, with a diffusion coefficient 
Din — o.^r (|27l . kll. |42|). For comparison with real systems, one-dimensional diffusion 
inside the channel must be matched to the three-dimensional diffusion outside the 



channel, through the choice of Tq (see e.g. (|27l . l28l l4ll . l55l)). For clarity, we re-derive 
this connection here without the inter-particle interactions inside the channel - see 
Fig. 6 for illustration. 

We denote the three-dimensional concentration of particles at the left side far away 
from the channel as c|f ; we assume that concentration on the right side far away from 
the channel is zero. At steady state, a density profile will be established such that 
the flux through the pore is F, the (three-dimensional) density at the pore entrance 
on the left is cl and the density at the exit on the right is cr] the corresponding 
one-dimensional densities are c\ = clPR"^, and c]^ — cjifiB? , where R is the channel 
radius, and /3 is a geometrical pre-factor that depends on the shape of the channel 
opening (/3 = tt for circular opening). 

At steady state, the flux that enters the channel from the left is (jHI): 



J = a{c^ - CL)RD,nt = F 



(19) 



where a is a geometrical pre-factor that depends on the shape of the channel opening; 
a = A for a circular opening (j4l[ ). Note that if all the impinging particles would go 
through the channel, the entering flux would be Jq = ac^RDout - the flux to a fully 
absorbing patch of radius i?(|4lh. However, even in the absence of jamming, not all 
particles go through - some of them return back, after hopping back and forth inside 
the channel, as reflected in the returned portion of the flux — acr.i?Z?r„,t (j41l ). 
The flux that exits the channel to the right is (|27l l4l[ ) : 



Jout = acRRD.nt = F (20) 
The flux inside the channel, for a flat potential proflle, is: 



F 



LZ 



'-D- 



(21) 



where Z = (e^) is the average inverse Boltzmann factor of the attractive energy inside 
the channel, E <Q. Solving the above equations, we get: 



F 



Jo 



9 I a _L Pout 7 
^ 13 R Din 



And thus the fraction of the transmitted flux is 

1 



Ptr 



9 I Q L Dput 7 



On the other hand, equation ([8]) gives without jamming (J — > 0) 

1 1 



2 + {N-l)ro/r 2+^^ 



(22) 



(23) 



(24) 
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Finally, choosing 7'o/r 



/An ~ 



fjy^j^Z, the discrete and continuous formu- 



lations become equivalent as long asiV = L/a^l ()27l . 

The distance between the 'sites' models the excluded volume interactions between 
the particles. In this paper we make the most parsimonious choice: the distance 
between sites is equal to the size of the particle. This choice adequately captures the 
essential properties of hindered diffusion in narrow channels 



5lD- In 



principle, in some systems the actual 'diffusion step' can be smaller than the particle 
size. However, in known cases, the results remain qualitatively the same after proper 
re-scaling of the transition rates (0, [64, 6^, 6^, 67). We also found that the quality of 
the fits in Fig. 5 and overall conclusions are not sensitive to small variations in the 
estimates of the parameters of the model (data not shown) . 



Expressions for M 1 

For completeness, we present here the expressions for the general case 1 < M < N/2, 
Um = 1 The optimal exit rate (for the values of J, Af and N when the optimum exists): 



^max , Jlr{l - M) V J/r(-2M + iV + 1) 
° ' J/r{M -If + {2M - N -1) ^ ' 



The channel occupancy 

J/r(2 + (iV - \)rolr){N + {M - l)(iV - M)ro/r) 

2(ro/r(2 + l)ro/r) + J/r(l + (M - l)ro/r)(l + (TV - M)ro/r)) 

and the saturation current in the J jr — > oo limit: 

ro(l + (Af-l)ro/r) 

(M - \)[N - M){ro/rY + (iV - l)ro/r + 1 



(26) 



(27) 
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Figure Legends 

Figure 1. 

Schematic diagram of transport through a channel A. Schematic illustration 
of the transport through a narrow channel. B. Kinetic diagram of a one-site channel. 

C. Kinetic diagram of a two-site channel. 

Figure 2. 

Kinetic diagrams of transport through a channel of an arbitrary length A. 

Symmetric channel consisting of N positions (sites) . The particles enter the channel at 
a site M with an average rate J. B. Equivalent energetic diagram in the case when the 
exit rates are determined by the interaction (binding energy) with the channel. The 
exit rates at the channel ends are given by Arrhenius-Boltzmann factors of the energy 
barriers at the exits, and E^: ~ exp(— i?_,/fcT) and ^ cxp{—E^/kT) C. 
Equivalent geometry of the channel in the case when the exit rates are due to spatial 
bottlenecks at the channel ends. 

Figure 3. 

Efficiency of transport through a channel of an arbitrary length A. Transport 
efficiency as a function of the exit rate for J/r = 0.01, Um = 1 for different entrance 
sites M. Black line: M ~ 1,N ~ 10, gray line: M ~ A, N ~ 40; corresponding dashed 
lines show the probability of a particle to traverse the channel; it is identical to a single 
particle transport efficiency in the limit J — > - cf. text. B. Transport efficiency as a 
function of channel length A^, for the optimal value of exit rate Tq = {Jr/(N — 1))^^^, 
M = 1, J/r = 0.01, rim = 1- Black line: J/r = 0.01, dashed line J/r = 0.1 C. 
Transmitted flux Jo^t/ J^^t ~ equations (8) and (16), as a function of the normalized 
incoming flux J/r; black line: ro/r ~ 0.01, dashed line ro/r = 1; M = l;n„i = 1. Note 
that the transmitted flux saturates to a constant value J^^j in the jammed regime. 

D. Optimal exit rate as a function of the channel length for M = 1, J/r — 0.01, 
n„i = 1, (black line). Dashed line: same for J/r = 0.1. 

Figure 4. 

Occupancy of the channel at the jamming transition A. Occupied fraction 
of the channel at the jamming transition, Tq = fmax, as a function of the channel 
length N, for different values of the incoming flux J/r. It shows that the channel can 
be occupied to a considerable degree - up to half of the available sites - before the 
jamming becomes significant. B. Densities at the entrance site 1 (black line) and exit 
site N (dashed line) as a function of the incoming flux J/r for ro/r = 0.1, = 5, 
rim = 1- Density at the entrance site saturates to 1, which causes the saturation of 
the transmitted flux. Density at the exit site stays low even in the regime when the 
transmitted flux through the pore saturates. 

Figure 5. 

Flux through nano-channels: comparison with experiment A. Flux through 
the nano-channel as a function of the outside concentration of the transported ssDNA. 
Black dots - experimental data from Ref.|l6h for a nano-channel without trapping 
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inside. Corresponding black line - theoretical fit from eq. ([8]) with Um — 6, Z=l, 
-Din/fout = 0.42, N ~ L/ (25). Red dots - experimental data from Ref.l|l6l) for a nano- 
channel with ssDNA hairpins grafted inside the channel, which are complementary to 
the transported ssDNA. Corresponding red line is the theoretical prediction of eq. 
dH) with n„ = 3, Din/Dout = 0.0042, Z = 0.00007, N = L/{2S) B. Reduction of 
the flux through the channel as a function of the number of mismatches between 
transported ssDNA and the ssDNA hairpins grafted inside, relative to the flux of the 
perfect complement ssDNA measured at the feed ssDNA concentration 9//M. Dots - 
experimental data from Ref. (flil ) for a single mismatch at the edge of the transported 
DNA segment; square - single mismatch in the middle of the transported ssDNA 
segment; line - theoretical model; same parameter values as used in panel A - cf. text. 

Figure 6. 

Three-dimensional diffusion outside the channel Schematic illustration of the 
three-dimensional diffusion outside the channel and one-dimensional diffusion inside. 
See text in Appendix. 
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A Channel is represented by N positions (sites) 
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